ITEP-TH-32/10 
MPP-2010-167 



Low-Energy Theorems from Holography 

Johanna Erdmenger^, Alexander Gorsky *, Petr N. Kopnin**, Alexander Krikun**, 

and Andrew V. Zayakin 1 '* 1 ' 1 

* Institute of Theoretical and Experimental Physics, 
\ B. Cheremushkinskaya ul. 25, 117259 Moscow, Russia 

' Max- Planck- Institut fur Physik (Werner-Heisenberg-Institut), 
■ Fdhringer Ring 6, D-80805 Miinchen, Germany 

* Moscow Institute of Physics and Technology, 
Institutsky per. 9, 14 1 700 Dolgoprudny, Russia 

Q-i| b Fakultdt fur Physik der Ludwig- Maximilians- Universitat Miinchen und 

^ , Maier- Leibniz- Laboratory, Am Coulombwall 1, 85748 Garching, Germany 

* Dipartimento di Fisica, Universitd di Perugia, I. N.F.N. Sezione di Perugia, 
^ \ Via Pascoli, 1-06123 Perugia, Italy 

°^ ; Abstract 

In the context of gauge/gravity duality, we verify two types of gauge theory low-energy 
theorems, the dilation Ward identities and the decoupling of heavy flavor. First, we 
provide an analytic proof of non-trivial dilation Ward identities for a theory holo- 
graphically dual to a background with gluon condensate (the self-dual Liu-Tseytlin 
background). In this way an important class of low-energy theorems for correlators 
of different operators with the trace of the energy-momentum tensor is established, 
which so far has been studied in field theory only. Another low-energy relationship, 
the so-called decoupling theorem, is numerically shown to hold universally in three 
holographic models involving both the quark and the gluon condensate. We show 
this by comparing the ratio of the quark and gluon condensates in three different 
examples of gravity backgrounds with non-trivial dilaton flow. As a by-product of 
our study, we also obtain gauge field condensate contributions to meson transport 
coefficients. 
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1 Introduction 

On the long road towards a holographic description of QCD, there are some milestones 
corresponding to exact relations which have to be satisfied also in any holographic 
model. These are so called low-energy theorems [1] (see e.g. [2] for review). In field 
theory these are statements which impose restrictions on the various correlators. The 
purpose of this work is to compare holography to field theory by considering the low- 
energy theorems concerning one- and two-point functions of a strongly coupled gauge 
theory on both sides of the correspondence. We report nice non-trivial agreement in 
two important cases: the dilation Ward identities and the decoupling theorem for the 
heavy flavor. Recently the validity of a related class of theorems (QCD sum rules) 
was shown holographically in [3] at finite temperature. Apart from demonstrating 
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the validity of low-energy theorems, a particular result of our analysis is a statement 
on the IR universality of theories dual to three scale-dependent backgrounds with 
non-trivial dilaton flow. 

First, we aim at realizing the QCD low-energy theorems explicitly, for instance 



in the conformal case: The right-hand side is expected to be zero in a conformal field 
theory where all condensates vanish. For an explicit expression for the correlators 
of energy-momentum components see e.g. [4]. Thus for a nontrivial test we need a 
background which is different from AdS in the IR, dual to a non-conformal field theory, 
for instance with a gluon condensate. There are a number of models which generalize 
the original AdS/CFT correspondence to the backgrounds corresponding to non- 
vacuum states of M = 4 SYM or to non-conformal and non-supersymmetric theories. 
We use the self-dual background by Liu and Tseytlin [5] with non-zero expectation 
value of the gluon operator (trG 2 ) in this part of our work. To perform the test 
of dilation Ward identities, we calculate the two-point correlators (trG 2 (x)trG 2 (0)), 



(trG 2 (x)trGG(0)), (T(x)trG 2 (0)), (T(x)trGG(O)), (T^(x)T^(O)) in this background. 



The analysis of correlators is easily performed for non-zero frequency. In this way 
we reproduce the results for transport coefficients, extending the analysis to the case 
of the non-conformal backgrounds considered. First of all, we calculate the rj/s ratio 
of shear viscosity over entropy via (T xy T xy ), which was performed for the conformal 
case in [6, 7, 8, 9]. Here we find using suitable holographic renormalization that 
condensate corrections to -|t->o = j- are absent in the Liu-Tseytlin background, i.e. 
the nonzero VEV of the gluon field strength (trG 2 ) does not affect the value of t]/s. 

Secondly, we check the relationship between two-point and one point functions 
in gauge theory with fundamental fermions, known as decoupling relation 



Fundamental fermions are introduced in our system via probe D7 branes, see e.g. [10]. 
The D7 branes represent the fundamental degrees of freedom, being convenient loca- 
tions for the fundamental strings to end and to be thus endowed with a global SU {Nf) 
flavor symmetry in the Maldacena limit. The length of the strings corresponds to the 
quark mass, and the subleading term in the asymptotics of the embedding coordinates 




(1) 



where T = is energy- momentum trace on the boundary. This is trivially satisfied 




12m(qq) . 



(2) 
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to the condensate. A non-trivial test of the theorem considered is possible only for 
an IR-non-trivial metric. For that purpose we use three different dilaton flow back- 
grounds with gluon condensate: the self-dual Liu-Tseytlin background mentioned 
above [5], the Gubser-Kehagias-Sfetsos background [11], [12] and the Constable- 
Myers background [13]. All of these are examples for non-trivial dilaton flows. A 
remarkable universality among the three models and agreement with standard field 
theory is observed. 

Let us now present the two low-energy theorems discussed in this paper. 



Dilation Ward Identity. It was argued in [1] that the following dilation Ward 
identity holds within field theory, 

limi J e tqx d A x (t ^^trG 2 (0)|^ = (-d)(0) [1 + mass-dependent terms] , 

(3) 

where d is the canonical dimension of the operator O, T{-, ■ } stands for the time 
ordered product and the one-loop beta-function is normalized as 0(a s ) = — b = 
yiV c — §iV/. Identities for higher correlators are also available: 

i 2 j d 4 xd 4 y(T^O(x),^^-tiG 2 (y),^^-tiG 2 (0)^ = (-d) 2 (0) [1 + mass-dep.] . 

(4) 

For the gluon field strength operators we obtain: 



/(r{^ W ,^,o,}>^(^ 



(5) 



Decoupling Theorem. Novikov, Shifman, Vainshtein and Zakharov derived in 
[1] the following equation for light quarks by considering the regularity of the beta 
function 

= -£<*>. (6) 



dm q \ 7r / b 
This low-energy theorem for heavy quarks is recovered also in an independent manner 
in [14]. Besides, for heavy quarks the following relation due to Shifman, Vainshtein 
and Zakharov holds 

m(gg) = -l(^trG 2 }. (7) 

The derivation of this equation is found in [15]. It expresses the continuity of the 
energy-momentum trace at the flavor number thresholds of the beta-function. The 
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factors 12 and 24 in the equations above are universal, they do not contain N c or 
Nf. In this paper, we shown that relation (7) holds holographically in the three 
dilaton-flow backgrounds to great accuracy. 

A related calculation, the holographic derivation of the Veneziano-Witten for- 
mula relating the mass of the rj meson and the topological susceptibility of pure 
Yang-Mills theory, was performed in [10]. The holographic conformal anomaly was 
previously considered under finite temperature in the 5-dimensional model with a 
dilaton potential adjusted in such way that both confinement and the correct UV 
behaviour of the coupling are reproduced [16]. 

This paper is organized as follows. In Section 2 we describe technicalities related 
to finding correlators. In Section 3 we describe the holographic description of the 
dilation Ward identities. Section 4 contains our main result - the derivation of the 
nonperurbative decoupling of the heavy flavor in the different dilaton flow models. In 
the last section we discuss the importance of having established the decoupling and 
scaling theorems holographically. Several necessary facts concerning the models are 
collected in Appendix A, while Appendix B concerns the derivation of the transport 
properties of the models under consideration. 



2 Recipes of AdS/CFT 



For later use, let us briefly review the AdS/CFT prescription for calculating two- 
point functions, emphasizing in particular the derivation of the gauge-fixing and the 
Gibbons-Hawking term. In the analysis of the boundary term we follow here very 
closely the analysis of [4]. A reader familiar with these technicalities can proceed 
directly to the next section. 

We consider the general rules for two-point functions and calculate the matrix 
of correlators 



Mij = (OiOj)]^ 



5 2 S 



full 



The standard wisdom on finding Green function of the fields present is to set the 
action of the type 

S b uik = J d 4 xdz<P' 2 g zz ^ (9) 
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out onto the boundary as 

Sboundary = I (Px^fi g"* y/g~\ z -> . (10) 



The correlator in terms of bulk-to-boundary Green functions G(x, z) of the field <fi is 
given by 

(O(x)O(0)) = G(x,z)d z G(0 : z)\ z=0 . (11) 

In our case two additional difficulties arise. First, the correct boundary term should 
be supplemented by the Gibbons-Hawking term [4], which makes a theory defined 
on manifold with boundary globally diffeomorphism-invariant. Second, the bilinear 
action of fields' fluctuations is non-diagonal, this means that we shall be dealing with 
a matrix of Green functions rather than with separately-treatable ones. 

Let us define Green function matrix. Namely, if field $j has a bulk solution <&i(z), 
satisfying z <5i $ i (z)| z _ 5 . = $j, then by definition 

Let us establish the correct boundary term. The full action of our bulk theory is 
actually [4] 

Sfuii = Siod + Sdiv + S i( i (13) 
where the Gibbons-Hawking term 

Sm = ~2d z J d A x^l - c J d*x^f 4 , (14) 

is here given by 

04 = det(gij), i = 0,1,2,3. (15) 

The constant c can be fixed arbitrarily to our convenience, e.g. as in eq. (4.15) in [4]. 
The other piece which one has to take into account is the full divergence term Sdi V , 
which does not affect equations of motion, but does change the appearance of the 
action and makes it diagonal in terms of physical degrees of freedom of the graviton. 
It is the well-known fluctuation term 

3 
2 

the vector W is (see [17], Vol.11, §96) 



Sdiv = ~d,W», (16) 



-9(g a "ST^-g^5T%), (17) 
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where 5T^ = T^(g + h) — r^(^). This constitutes the gauge-fixing prescription for 
our problem. 

Consider now the second variation of these actions in fluctuation fields; denote these 
second-order expressions as S^ d , S d f v , S^J respectively; they contain both fields and 
their derivatives. The two-point correlator is then 



d 2 C 

(0,0,) = K ik d z K jm + K ik 

k m 



f)2 c( 2 ) f)2 c( 2 ) 

d<s> k d$ m K ^ 



d<$> k d$ m 



here L is Lagrangian density of the bulk action: 

Sbuik = S[J d + S div = I dz C 



(19) 



The above structure is obvious from the following reasons. Consider the bulk action 

,20 5$ m (z) 5 2 S bulk 5$ k (z) 

Obuik — — — — r= , {ZV) 

where 



5 2 S 



bulk 



dz 



8 2 L 



7 d z 6$ m dz5®k + 



d 2 L 



d$ m dV k 



5$ m d z 5<f> k + 



d 2 L 



d® m d<$> k 



-6<$> m 6$ k 



(21) 



Taking into account that Green functions of field fluctuations by definition satisfy 
equations: 



a d 2 L d 2 L 



d 2 L 



5$,(z) 



0. 



(22) 



one sees that the only contribution of Sbuik into the correlator will be, after taking off 
the derivative and integration, the term: 



d 2 L 

5 2 S bulk = 8$ m {z) d z 5$ k (z). 

m k 

Now remembering the definition of Green function matrix 



5$, 



(23) 



(24) 



we arrive exactly at (18). Then there is the purely boundary term (Hawking-Gibbons 
term). It does not require the above procedure, since it already sits on 4d. Then it 
contributes the following: 

S 2 S M = d,6$ m 6$ k + _f ^ 6Q m 6* k . (25) 

The action contains no more than one derivative term, which is due to normal 
differentiating of extrinsic curvature, thus J|t| = 0. This contributes the other two 
terms into the correlator (18). 



3 Low-Energy Theorems 

In this Section we calculate the matrix of the two-point correlators for the gluonic 
operators and components of the energy-momentum tensor. Then we compare these 
to one-point correlators and find that the correct scaling relations from field theory 
are satisfied on the gravity side. We begin by introducing the Liu-Tseytlin model in 
which we will perform our calculations in this section. 

Liu— Tseytlin model. In the Einstein frame the bulk action of IIB superstring 
theory is [5] 

Sio = ^y^l J dW ^o (R ~ \{dAf - \**{%Cf - i|F 5 | 2 ) , (26) 

where R is the curvature, is dilaton, F 5 is 5-form and C is axion. 

The Liu-Tseytlin model is a generalized background for holography those which 
possesses self-duality. It describes a field-theory flow from a strongly-coupled con- 
formal theory in the UV to a theory with condensate trG 2 in the IR. By virtue of 
self-duality it is still supersymmetric. However, it possesses a scale parameter, which 
makes it closer to real-world physics. The self-duality is provided by the presence 
of a non-trivial axion field. Despite the presence of the scale, it is conformal in the 
UV; in the IR the dilaton singularity is determined by the gluon condensate trG 2 . 
Within supergravity this background is understood as "smeared" D(— 1) brane with 
a usual stack of D3-branes. Since D(—l) brane is an instanton in 10D, the resulting 
4d theory can be considered as having an instanton-gas type of vacuum, which is ad- 
vantageous for QCD purposes. Moreover, this background is confining (in the sense 
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of Wilson loop linear behavior at large temporal separation), and the string tension 
is proportional to the condensate. Of course, we do not claim to produce any real 
QCD results in this framework, but we believe it to be a very useful toy model. 

For the Liu-Tseytlin background [5] metric in Einstein frame looks like the 
standard conformal solution 

ds 2 = g^ v dx^dx v = R + y/h 3 M , (27) 

but the dilaton is modified by the smeared instanton (nonzero density of D(—l)) 

= h- lt (28) 

and an axion is present 

C = J- - 1; (29) 
the D3 and D(—l) form-factors are: 

h = z\ (30) 

and 

h^ = l + qz A . (31) 

The parameter q is the crucial quantity for us, since it measures the degree of IR- 
non-conformality of the theory (remember that in the UV, the theory is conformal 
and its /3- function is zero). 

The Tseytlin-Liu background has been successfully used for a number of appli- 
cations, e.g. calculating meson spectra [18, 19, 20, 21]. In all these applications, its 
relevance to QCD has been demonstrated. In [22] a finite-temperature extension of 
the [5] solution has been found, which has been a further motivation to apply it to 
realistic high-energy quark-gluon plasmas. We shall employ Liu-Tseytlin background 
to test dilation Ward identities in Section 1 and decoupling relation in Section 4. 

Holographic normalization of the operators 

Here we consider normalization of the gluon field strength operator; the normal- 
ization of the quark operators will be considered in the next Section. According to 
the AdS/CFT dictionary we state that the fluctuation S(f>(z,Q) of dilaton field 

<j>(z,Q) = M*) + fy(z,Q) (32) 
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is dual to the operator 0$, proportional to the QCD scalar gluonic operator 

tr (G 2 ) = lo . (33) 
We can fix the normalization constant c$ by comparing the two-point functions 

(0,0,) = cJ(tr(G 2 )tr(G 2 )). (34) 
At large momenta the leading behavior of gluonic correlator in QCD is [23]: 

N 2 - 1 

<tr(G 2 )(Q)tr(G 2 )(Q)> = -^^Q 4 ln(Q 2 e 2 ). (35) 

To obtain a two-point function from holography we take the second variation of the 
action computed on a classical solution. In the vicinity of the boundary of AdS§ the 
action (26) for the fluctuation is: 



^ = R l ,4 / d 4 xdz\\ [-{dMf ~ d^8<j> + 2e^8<j>{d z C) 2 ] . (36) 
g^ZTTj'a^ J z 6 2 

Here we have taken the near boundary limit r ^> L (so that r 2 ~ p 2 ) and changed 
coordinates z = % 3 is the volume of the S5 sphere, R 8 came from the determinant 
of the metric (y/g = ^-)- The last term containing the profile of axion field is 
negligible at the boundary (small z) because d z C(z) ~ z 3 . We can find the bulk-to- 
boundary propagator of <p(z, Q) at small z and large Q 2 . It is 

tp(z, Q) = ^K 2 (Qz), <p(0, Q) = 1, (37) 

where Ki is McDonald function of the second kind. Now we can compute the second 
variation of the action. It is 



5 2 S d n 3 R 8 1 . _,0^(z,Q) 



Wo £ 2 (2tt)V 4 2^ v z 



3 



N 2 1 



Q 4 ln(Q 2 e 2 ), (38) 



=e 4(2tt) 2 8 

where we used the definition R A = 47rg s a' 2 N c and the asymptotic of McDonald func- 
tion. Comparing this result with the expression of QCD we find 

= ^tr(G 2 ). (39) 
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To establish a relation between gluon condensate and the expansion coefficient of 
the dilaton field we compute the vacuum expectation value of 0$ at zero momentum 
taking the first variation of the action with respect to the boundary value of the field 
O . At zero momentum near the boundary the dilaton field behaves as 

<f>{z) = O + faz 1 . (40) 
For the dual operator given by (39) we find 



SS d tt 3 R 8 1 , n M{z,Q) 



^ 2 (27r)V 4 2 rv z 3 



N 2 

c 



4(2tt) 



(41) 



From (39) and (40) we get the expression for the gluon condensate 



(tr(G 2 )> = 4^20, = iV c 2 ^04. (42) 

In the Liu-Tseytlin model the infinitesimal fluctuations of the fields on the bulk 
couple to the operators trG 2 , trGG, T^ u in the boundary M = 4 SYM theory. More- 
over, in the Liu-Tseytlin model the dilaton field behaves as = 1 + qz A , so the 
parameter of solution 04 in (42) equals q and the scalar and pseudoscalar gluon con- 
densates are nontrivial and equal to the value given in (42), i.e. 

(trG 2 ) = (trGG) = tfj^q. (43) 

Correlators at Zero Frequency Fluctuation terms are defined as 

= 4>c + <p, 

C = G + £, (44) 
9 = 90fiv + h/iv 

We consider the following interaction term to provide a correspondence with the 
boundary theory: 



Sin.t — I d X 



- e 



(45) 



which, after introduction of useful self-dual and anti-self-dual components 



G ± = (46) 
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and splitting axion and dilaton fluctuations into a new couple of variables 

ri k = tp±Z, (47) 

becomes 



Sint. — I d X 



l -T, u h^ - tjL (fj + txG+ 2 + fj- trCT 2 ) 



(4c 



Here bars denote four-dimensional sources, which are proportional to boundary values 
of five- dimensional fields: 

hfj, u = z 2 h^ u \ z=Q , ff = ?7 ± | z= o, (p = <f\z=o- (49) 

Fluctuations of F5 are fully determined by thus there is no independent source 
for them. 

Let us choose the gauge h 5 ^ = 0, k^h^ = 0, u^h^ v = 0, where wave-vector 
k = (u, 0, 0, k), constant vector u is u = (1, 0, 0, 0). We work with five fields: 

$i = (?7 + , h u + h 22 , hi - h 22 , hi2, Tf), (50) 

z = 1, ... 5, each coupled to the corresponding Oi operator 2 

ArG+ 2 1 3 1 1 1 trG" 2 \ 

Oi - [~^, ~T, , gT n - -T 22 - -T33 - -Too, ^, -^=- j , (51) 

with G + and G~ the self-dual and anti-self-dual parts of G, respectively, via 

5 



S in t= / d A xdzy^ j Oj^j. 
•* i=i 



(52) 



(53) 



The relevant part of the fluctuation action in the bulk is 

^double deriv. = J ^ + + ,^ + | # 

One should not be mislead by its diagonal structure; besides the diagonal terms with 
double derivatives, the full bilinear action contains terms which make it non-diagonal. 

The boundary Gibbons- Hawking action term is 

^.derivatives = J ^1 (^^2 + J^/^)/^) + ( C $ 2 (z) + 4^ a (z))) . 

(54) 

2 Some of these operators, e.g. the O3 are not of immediate interest; however, it costs no additional 
effort to incorporate them into the calculation, so we work the correlators out for them as well. 
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The full system of equations upon Green functions (22) in the given background (27)- 
(31) is cumbersome and therefore is given in the Appendix B, eq.(119). Note that for 
the F 5 form we always have 5F = — 2/r 3 $ 2; which solves automatically the equations 
of motion for this field and at the same time retains the constancy of the Ramond- 
Ramond flow J s5 F 5 = N c . 

It is instructive to start with zero-frequency correlators (setting u = in (119) 
in Appendix B). Subsequently, we introduce finite frequencies u. In this case we find 
oscillatory solutions (Bessel functions) (121) instead of the rational ones (120). The 
limit oj — > of the finite frequency result coincides with our previously found result 
at u = and thus provides an additional check of the validity for our procedure. 

The solutions (121) contain ten modes labelled by coefficients C,, i = 1 ... 10. 
One would expect that out of the ten modes five must be IR finite, yet quite unex- 
pectedly six there are six IR finite modes (Ci, C4, C5, Ce, CV, C9), and the remaining 
four are infinite. An extra constraint is therefore necessary to make the Green func- 
tion matrix (12) a well-defined 5x5 matrix. We require that the resulting corelator 
matrix be symmetric, which is equivalent to the condition C5 = Cq/2, which removes 
exactly one redundant degree of freedom. 

The Green function matrix is then (recall that = -^m)'- 





C(j> trG+ 2 


81" 


o 3 


04 


c trG- 2 


C<t> tiG +2 


qz 4 — qe 4 + 1 














8 A« 





1 











o 3 








1 

2 2 



















1 





c<j> trG~ 2 


-2q (e 4 - z 4 ) 











4 4 i 1 
qz — qe + 1 



(55) 



with the Oi as given by (51). q is the non-conformality parameter defined in (31). 

As a result, combining our knowledge of Green function matrix (55), the bound- 
ary action (54) and the derivative piece of the bulk action (53) we obtain the matrix: 



M 





trG +2 


iTt 1 

8 /■« 


o 3 


04 


trG 2 


c<p tiG +2 


-Aq 


-2q 








-2q 


8 C 


-2q 


1 

4e 4 











o 3 








1 

4e 4 



















1 

e 4 





trG~ 2 


-2q 















(56) 
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which contains information on the correlators of Oi, Oj via the following relation 



N 2 



16tt 



2 ±y±% r 



(57) 



Some comments are due here. The singular terms \ are expected due to the 
divergencies on the field theory side; they are subtracted by a holographic renormal- 
ization procedure, analogously to field-theoretical subtraction. The asymmetry in 
0\ -H- 05 is also expected: what we consider is a self-dual configuration, therefore, 
the self-dual and the anti-self-dual operators have different properties. 

Using the matrix elements obtained above, we can now establish the low-energy 
theorems. After normalization according to (39) we have 



d 4 x(tiG +2 (x)T(0)) = 4<trG +2 (0)) 



d 4 x(trG~ 2 (x)T(0)) 
d A x (trG 2 (x)trG 2 (0)) 
d A x (tiGG{x)tiGG{0) 



(5J 



0. 



where T = T vu . Here we see that the first and the second lines of the equations 
above (58) constitute exactly the statement of the low-energy theorems 

(OT) = dirn(0)(6). (59) 

Note that (tr {G~) 2 ) = 0. 

The third line of (58) must be compared to the field-theoretical result 



(trG 2 trG 2 ) ~ l/A^trG 2 ), (60) 



where in standard perturbation theory, j3 Q is the one-loop coefficient of the beta- 
function. This equation reflects a breaking of the conformal symmetry. For the 
Liu-Tseytlin model the standard beta function vanishes. Nevertheless, the massive 
parameter q generates additional terms in the effective action. This gives rise to the 
contribution T^a ~ trG -2 to the trace of the energy-momentum tensor at the operator 
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level. This is consistent with the low-energy theorem given by the third line of (58). 
On the other hand, (trG~ 2 ) = 0, thus the expectation value of the energy-momentum 
tensor and the vacuum energy vanish, ensuring consistency with supersymmetry. 

The fourth relation in (58) implies that the topological susceptibility of the 
vacuum, which is proportional to this correlator [24], vanishes in the Liu-Tseytlin 
model, which is in the agreement with the fact that the model is supersymmetric 3 . 



Correlators at Finite Frequency Now let us analyze the finite-frequency solu- 
tions. The solutions are given in Appendix, eq. (121); only relevant modes shown. 
Unlike the u = solutions, which were exact solutions, here <&2{z) and $5(2) are 
powerlog expansions in u and r. Since we are interested in the near-UV behaviour 
of Green functions, and eventually expand correlator matrix in powers of u, this 
approximation is reasonable. The matrix of correlators becomes: 



M 





C, trG+ 2 


177* 

8 C 


o 3 


04 


c> trG- 2 


04, trG+ 2 


-Aq 


-2q 








log(we)w 4 r, 


8 A* 


-2q 


log(we)a; 4 
32 



















log(we)oj 4 
32 



















log(a;e)cij 4 
32 





cp tiG~ 2 


log(we)w 4 ^ 















(61) 

The most interesting physical implication of this correlator matrix comes from the 
(T xy T xy ) element. It is proportional to -\t=o, and here we observe its independence of 
q. This fact is not trivial from dimensional considerations, since we possess another 
dimensionful parameter, the frequency u. Thus we have established 



V 



T=0 



1 

471"' 



(62) 



As a bonus of this calculation, in the Appendix A we easily elaborate the matrix of 
quarkonium transport coefficient based on the above correlator matrix. 



3 Note that in the D4/D6 model [10] the topological susceptibility does not vanish. However there 
is no contradiction between these facts, since the model of [10] breaks supersymmetry (similarly to 
Sakai-Sugimoto model), whereas Liu-Tseytlin model retains supersymmetry. 
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4 Holographic Decoupling of the Heavy Flavor 



4.1 Physics of Decoupling 

In this Section we holographically derive the central result of this paper, which is 
known as "decoupling relation". In can be found in [15]: 

^(G%G%) = -12m q (qq). (63) 

The derivation of this relation is somewhat intuitive, but let us still restate the argu- 
ments by Shifman, Vainshtein and Zakharov. For vacuum expectation values of the 
different operators pertinent to light quarks the parameter of expansion is quark mass. 
For heavy quarks we expand in the inverse quark mass and set external momentum 
to Q 2 ~ 0. Let us suppose there exists a quark for which both expansions, small and 
large m are true. As it is in particular a "heavy" quark, the quark condensate can 
be done perturbatively from the triangle diagram with gluons as "vacuum sources" , 
shown in Fig. (1). 




Figure 1: Vacuum diagram with heavy quarks depicting (qq) as gluon-driven quantity. 



One can understand the argument from which the relation (7) emerges as follows. 
Consider the trace of energy-momentum tensor of a gauge theory. For low quark 
mass there is beta-function contribution from the quark, for heavy quark there is only 
the gluonic contribution to the beta-function, yet there is quark chiral condensate is 
present: 



— Nr~- 



3 J 8tt 



above threshold, 



(64) 



11 



a 



— N c — trG + mqq, below threshold. 



>7T 
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When the two are equated at some intermediate scale, the necessary relation (7) ap- 
pears. Equating small and large m domains happens on the ground that we select the 
scale at which the heavy quarks "decouple" from the one-loop polarization operator. 
Hence this theorem is also known as decoupling relation. A picture of condensate as 
function of quark mass is given in [1]. 



4.2 Decoupling in Specific Backgrounds 

We now establish relation (63) holographically by considering different backgrounds, 
those of Constable and Myers [13], of Gubser [11] and of Liu and Tseytlin. The Liu 
and Tseytlin background (27) was already discussed above in the Section 1. The 
Constable — Myers background in the Einstein frame has the metric 



i 



ds 2 = ± -=L dxl + ^ (dr 2 + r 2 dtt 2 ) 

rfr A* -4 \ 5J 



( 6 4 +r 4 W ( r A h A\ ( b 4 +r 4 N > 4 1 2 W*)/ (b*+r~ 



where 



(65) 

(^4 (66) 



and the dilaton is 



ft 4 _|_ r 4\ 5V 10 "3^ 



H - V.' ' (67) 

axion is zero, and F 5 = £5^, where €5 is the unitary antisymmetric tensor in the S5 
directions. 

The chiral condensate and meson spectrum involving a Goldstone boson were 
obtained in [25] by embedding a D7 brane probe into a Constable-Myers background. 
Masses of heavy-light mesons in this background in D7 model were obtained in [26]. 
The quark condensate, pion decay constant and the higher order Gasser- Leutwyler 
coefficients were calculated for D7 model in this background in [27]. D7 embeddings 
were argued to be stable in this background [28, 29]. 

One of the first non-conformal backgrounds introduced into AdS / CFT was con- 
sidered by Gubser [11]: 



ds 2 =\ll-^r 2 dxl + ^(dr 2 + r 2 dnl), (68) 



16 



dilaton in this background is 




and the axion is zero. Originally it was intended to model confinement, yet it 
became also useful for introducing the gluon condensate. Shortly before Gubser, 
this background was also obtained by Kehagias and Sfetsos [12] in a less convenient 
parametrization. 

Introduction of fundamental fields. We are modelling the fundamental fermionic 
degrees of freedom by embedding the D7 brane into one of the three backgrounds de- 
scribed above. The Dirac-Born-Infeld action for the D7 brane embedding in Einstein 
frame is given by 

Sm = I* (70) 

The embedding of D7 is made as shown in the following table: 



AdS 5 x S 5 





1 


2 


3 


4 


5 


6 


7 


8 


9 


D7 


+ 


+ 


+ 


+ 


+ 


+ 


+ 


+ 







(71) 



One can get an image of the corresponding physics in Fig. (4.2), where string modes 
generating specific sectors of the spectrum are shown. 




3-7 chiral multiplet Q 
7-3 chiral multiplet Q 
3-3 N=4 vector multiplet 



Figure 2: Scheme of the D3-D7 geometry and the corresponding string/field modes. 

We look for embeddings of the form 

X 9 = w(p), X 8 = 0, (72) 
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where embedding function w, worldsheet coordinates and target space coordinates 
r, p are related as follows 



w 2 {p) = r 2 - p 2 , 

p = + a + e 7 + & 1 ] 

Nf quark flavours can be considered introducing Nj corresponding D7 branes with 
embedding coordinates Wi,i = l...Nf. If the quark masses are equal, D7 branes 
form a stack and the action (70) is multiplied by the factor Nf. In the following we 
restrict ourselves to the case of just one flavour for simplicity, considering only one 
embedding coordinate w(p). Using these definitions we easily construct the equations 
of motion for w(p), 

2pg m {r)w'{p) (w'(p) 2 + 1) <4(r) - 2w(p) (w'(p) 2 + 1) (g 55 (r)g' Q0 (r) + goo(r)g' 55 (r)) + 
+£550) (2pg' 00 (r)w'(p) 3 + 2pg' m {r)w\p) + rg 00 (r)w"(p)) = 0, 

(74) 

where the corresponding ga should be taken for each respective metric. We solve them 
numerically at different values of the vacuum parameters and fields, corresponding to 
the boundary conditions at p — > oo; a typical embedding is shown in Fig. (3). 




fa] Constable-Myers (bj Gubser-Kehagias-Sfetsos 



Figure 3: Typical embeddings of D7 branes. 



4.3 Normalization of the "Quark" Operators 

Following the same steps as in Section 3 we explore the scalar field w dual to the 
operator qq, where q is the quark field. It is described by the action of the D7 brane 
(70), for which Wi is embedding coordinate. Here and after we are dealing only with 
flavour i and will omit this index where it is possible. The action for the fluctuations 
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of w is 



S R 



2ti 2 R 4 



<7 s (2tt)V 4 



d 4 xdze* 



2z 



(d z w) 



— d a wd^w 



(75) 



Here we change coordinates the same way as in (36), 2ir 2 is a volume of 3-sphere R 4 
comes again from the determinant of the metric \fg® = nr- In the limit of large 
momenta near the boundary the bulk-to-boundary propagator is 

w(z, Q) = Qz Kx{Qz), w(0, Q) = 1. (76) 

The scalar field is dual to the operator O w , which is proportional to qq = —O w . We 
compute two-point function of O w to fix the normalization 



{O w O w ) 



S 2 S 8cl 2n 2 R 4 A ^.d z w(z,Q) 
owqOWq g s {27r)'a' q: 2 z 



(77) 



2 2\ 



-Q'HQ'e 



2(2vr) V 2 2 

Here the fact is used that ^boundary = 1 [5] and again R A = Ang s a' 2 N c . We compare 
this result with the QCD calculation (see eq. 4.27 in [15]), 



(qq qq) 



167T 



2,2\ 



:Q 2 ln(Q 2 e 



and find 



i 



2ira 



- f qq- 



(7* 



(79) 



At this stage we can identify the boundary value of the field wq = w\ z= q. It is the 
source of O w = c w (qq), so it is proportional to the quark mass wo = —M. Thus we 
have 



M 



2na 



(80) 



To identify quark condensate (qq) we compute the expectation value of O w at Q = 0. 
In this limit near the boundary the supergravity field takes the asymptotic form 



w(z) = w + w 2 z 2 . 



The result is 



owq g s {2rc)'a' i 2 z 



2(2vr) 4 a' 2 



2wo. 



(821 
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The quark condensate is normalized as follows 



^ = b ^ = w^ w * (83) 



To check the decoupling theorem, we have to study the relation 3 M ^ — for one 

specific quark flavour, i.e. Nf = 1. Expressed via the parameters of our model 
(42), (80), (83), it turns out to be 

M(qq) w w 2 



^(tr(G 2 )) N c 4V2a' 2 g$ M 04 



4tt 2 

with coefficients wq,W2,4>4, defined in (40), (81). It is convenient to express all coef- 
ficients via the expansion parameters in the coordinate r = — . We denote them by 
4> = 0o + pL w — a +^2- Obviously, these are related to the former defined in (40), (81) 
by </) 4 = jp,0J2 = Jr. Recalling that R A = <iirg s a' 2 N c and g\ M = <iiig s , we obtain 

M(qq) _ 1 ac ^ 



£*<tr(G*)> 4^2 V 
M{qq) 1 



4tt 2 

For the theorem 

M lnn\ 1 

(86) 



4^ 

4tt 2 



^(tr(G 2 )) 12 



to hold, the parameters a, 64, c must satisfy 

ac y/2 



bl 3 • < 87 > 

This relation, equivalent to the decoupling theorem, will be tested numerically below. 



4.4 Numerics 

We obtain numerically the dimensionless ratio of the solution coefficients |j, linearly 
related to the condensate ratio by (85). The ratio of the condensates obtained numer- 
ically is shown in Fig. (4) for the Gubser background. Similar pictures are obtained 
for the other two backgrounds. Each point in the parameter space represents an 
individual "measurement", that is, a solution for a D7-brane embedding at given 
gluon condensate and quark mass, from which the value for the quark condensate 
follows. By fitting the "experimental" points we estimate the value of the ratio and 
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a statistical error margin thereof. The mass is assumed to be the largest scale under 
consideration. 



tr G7 (m < \\t \\f >) 




We obtain numerically the following results for the dimensionless ratio of the 
condensates we are looking for: 



M(qq) 



Constable-Myers 12.0078 ± 0.005 
Gubser 12.25 ± 0.01 

Liu-Tseytlin 11.9192 ± 0.0020 



2 

Comparing the results to the correct analytic value — 4n2 M ^ q ) — = 12, we see agree- 
ment with good accuracy. The obvious universality of the three different metrics might 
signal that the decoupling theorem is insensitive to the details of the IR physics. 



5 Discussion 

Let us restate the main results of this study: 

• We have established a universal constant value for the ratio -^py in the holo- 
graphic duality with a good precision (0.5%), thus supporting the validity of 
the heavy flavor decoupling in holographic models of QCD. 
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• We have obtained a version of the low-energy theorem J (TO) = dim(0)(0) 
satisfied in holography with condensates for the pure glue sector. 

In addition we also find the following results 

• A non-trivial relation between two-point and one-point functions J (G 2 G 2 ) = 
const (G 2 ) has been established. 

• Shear and bulk viscosities have been shown to be independent of condensates. 

• The quarkonium diffusion coefficient has been obtained both at non-zero tem- 
perature and for a non-vanishing condensate in the Appendix A. 

The significance of establishing the decoupling ratio is its relevance to justifying 
phenomenological approaches to QCD via holography. We hope to find an analytic 
explanation of the amazing agreement which appears not to be a coincidence. Here 
we provide demonstration for a very important example of a statement which relates 
the quark and gluon sectors. This should encourage further development of realistic 
AdS/QCD constructions based on geometries with broken scale-invariance. 

Acknowledgements 

A.Z. is extremely grateful to Prof. Dr. Dietrich Habs for his generosity, encourage- 
ment and a fruitful atmosphere of Theory vs. Experiment interaction he has created 
in Garching providing a powerful inspiration for this work. Thanks to the Organizers 
of the International School on Strings and Fundamental Physics in Munich/Garching. 
A.Z. thanks Derek Teaney for correspondence on transport coefficients. The work was 
also supported by the DFG Cluster of Excellence MAP (Munich Centre of Advanced 
Photonics) (A.Z.), and by the Cluster of Excellence "Origin and Structure of the Uni- 
verse" (J.E.). The work of A.G. is supported in part by the grants PICS- 07-0292165, 
RFBR-09-02-00308 and CRDF - RUP2-2961-MO-09. The work of PK. is supported 
in part by the grant RFBR-09-02-00308 and by the Dynasty foundation. The work 
of A.K. is supported in part by RFBR grant no. 10-02-01483 and by the Dynasty 
foundation. The work of A.Z. is supported in part by the RFBR grant 10-01-00836 
supported by Ministry of Education and Science of the Russian Federation under the 
contract 14.740.11.0081. 



22 



A Quarkonium Transport in Self-Dual Background 



A.l Self-Dual Background at Zero Temperature 

Here we review the method of [30] for calculating quarkonium transport properties. 
The basic result of this discussion is a decoupled structure, in which the contributions 
of the fermionic part of the action will be separated from those of the gluonic part 
according to the pattern 



meson kin. coeff. 



two-point correlator 
(D3 contribution) 



(89) 



"meson mass shift 
(D7 contribution) 

Consider a complex field ip of a slowly moving meson of velocity v, coupled to some 
operators of gluonic sector, 



L = (f + vd t (f + ^2 CnV + O n (p, (90) 

n 

where coefficients c n are defined e.g. from D7 action of a dual model, which secures 
existence of mesons. The latter are understood as eigenmodes of fluctuations above 
the classical solution of D7 equations of motion. Interaction terms modify the spectra 
of eigenmodes in bulk; in terms of the boundary theory this amounts to meson mass 
shift. Coefficients C n are then introduced as "susceptibility" of mass with regard to 
switching on the operator O n : 

6M = -c n (O n ). (91) 
Considering one-particle dynamics we can obtain from (90) 

* = (92 > 

where 

Fi = / d 3 xtp+ViC n O n <p, (93) 



while correlator of two forces is directly related to transport coefficient 

K = - I dt{jr(t)F{0)}. (94) 



3 _ 

One can integrate field tp out of these relations and obtain finally 



K 



If 2T 

_ J k 2 d 3 kc 2 n —lm(O n O n )\ k , (95) 
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where 

(O n O n )\ k = J d*xO{t)j^ % *(O n {x)O n {G)). (96) 

Here the contributions of flavor dynamics and pure gluodynamics are decoupled; 
below we proceed in calculating the gluodynamical part (the two-point correlator); 
the coefficients c n being responsible for mass shifts are known in literature. 



A. 2 Self-Dual Background at Finite Temperature 

It is possible to obtain quarkonium diffusion and relaxation coefficients at finite tem- 
perature and condensate, extending the work [30] 4 to the background of [22]. This 
background has the metric 

,2 n2 /l-rVr 4 j2 dxl dr 2 \ n2 - 

ds = R { -2 dt " + 7# + ri(1 _^ ) J + R2dQ l (97) 

dilaton is 

e * = 1 + ^ l0g (r37W)- < 98) 

axion is related to dilaton in the same way as in the zero-temperature Liu-Tseytlin 
background 

C = e-* - 1. (99) 

Quarkonium transport coefficients are quantities which feel both the fermionic piece 
of the action (some embedded brane) and the gluodynamics. From the former comes 
mass susceptibility to condensate, from the latter - correlators of interest. In prin- 
ciple, it would make a good sense to work in a back-reacted metric, however this 
we postpone till the method is fully technically developed for the well-controllable 
Ghoroku-Liu-Tseytlin metric. 

For convenience we further use the variable 

u = rVT 2 , (100) 

which lives in the interval (0,1). We consider a reduced sector of the fluctuations, 
namely, those of fields rj~, hu+h22- The equations of motion are given in Appendix 
A (122). 

4 We thank Derek Teaney for providing us with his unpublished Notes. 
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We see now that the problem of fields coupling to each other is additionally 
burdened by presence of finite temperature. Yet diagonalization of these equations is 
possible by means of the following functional transformation 



7] + {u) 

h{u) 
fj-(u) 



7] + {u) 

h{u) + q(Cx- TV 2 log (1 - u 2 )) V + {u) 



qh{u) Fx 



l-tr 



V («)- 



(101) 



+q ( igbg 2 (i - it 2 ) - ^r 2) + C 2 ) V +(u) 



Now for each of the variables we can write down an equation similar to that for the 
simple dilaton modes: 



„. s u (u 3 + 6u + 4uj 2 + 4k 2 (u 2 - 1)) -3 , , 
V» + — , o, o 7T2 V( u ) 



Au 2 [W 



0. 



(102) 



,2 _ iy 

for which transport coefficient is known; we calculated it independently, and found it 



to be in agreement with the previous results [30] 



2uj 



TV 



(103) 



where C 7 = 4y (|) m 2.62. Knowledge of diagonalization matrix allows us to 
transform these results (at q = 0) into non- zero-condensate background: 



qA)- 



(104) 



where zero-condensate solutions are rotated to non-zero-condensate by the following 
rotation matrix in mode space: 

/ \ 



A = 7T 2 

\ 1/2/tt 2 J 
and the non-perturbed matrix of finite-temperature correlators is diagonal 



(105) 



(106) 



/ (trG +2 trG+ 2 ) \ 

(TT) 

\ (trG- 2 trG- 2 ) J 

whence one easily gets the mesonic transport coefficient by use of the following for- 
mula: 

E2 1 7r 

O 



2 d 3 k 2u 



32 J (2tt) 3 T 



(107) 
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where the respective mass susceptibility coefficients are obtained from considering the 
fermionic fluctuations coming from the embedded D7 brane piece of the action, and 
are defined via 

5M = -c {0), (108) 

where M refers to the mass of quarkonium. 

The correlators themselves are obtained in the following way, which we illustrate 
on the example of dilaton. We consider three domains: UV, IR and the intermedi- 
ate domain (we denote the latter QC for semiclassics, since semiclassical approximate 
solutions will be valid therein). The physical limitations are infalling boundary condi- 
tion on the horizon and reflected wave in the UV, which reduces number of unknown 
coefficients from 6 to 4. Then, we have matching conditions separate for each of the 
modes in the matching regions between UV and QC, an between QC and IR. This 
provides additional 4 constraints, thus the system is fully defined. In the UV the 
general solution to EOM is 



* = 2Uh ( 2 ^^ l + 2u (f - <J) K 2 (2V~uVB^) C, (109) 

Taking the UV asymptotic (tt — > 0) of 0, we see that physical boundary conditions 
are C\ — B,C% — 1, where B is related to correlator straightforwardly: 

= J™* (110) 

On the contrary, expanding it for large k, we get the form appropriate for matching 
with QC: 

d-2AV« 1.-9/2 .-5/4 

= e -2*V=^Ar 29 / 2 u- 5 / 4 - — • (HI) 



1-K 

The semiclassical equation has the approximate potential 



V QC = 7T~~~ 2~y ( 112 ) 
u (1 — u ) 

which allows to obtain the wave-functions in the standard way 

±Jpdx 

= -^p-, (113) 

where 

P = VVqc ~ E. (114) 
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The semiclassical solution near u = and u = 1 is 

' 2k ^ (e 4k ^A l + A 2 ) 



w-5-0 = 



(115) 

- 2m 

The IR solution with infalling boundary condition has only one degree of freedom: 



u — ^ i 



>IR 



f e ^c S c( TO + e -^ fc ^) . (116) 

Equating the QC solution branches with those of IR and UV solutions, we get 

ImB = ix 2 k 4 e- 2C ^ T } (117) 

as already stated above. Taking the integral over phase space (107) and performing 
linear transformation of correlator matrix (104), we get for transport coefficient 

i 6orm 6 

K = -T 9 [ct rG+2 (l + 2qn 2 ) + c T (l + qir 2 ) + c trG - 2 ] , (118) 

6 ^(i) 



where q are found in [30], c trG 2 = ^ (||) , c T = £ (jg^ , M being the 



meson 



mass. 



B Equations of Motion 

Here we shown the equations of motion for Liu-Tseytlin background in the graviton, 
axion and dilaton sector, corresponding to the pure glue sector on the boundary. The 
definitions of the fields are contained in eqs. (51)- (53). 
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' z ((quz 4 + u) 2 -32q 2 z 6 ]r)+(z)+ 

+ (qz 4 + 1) ((llgz 4 + 3) rf + (z) - z (qz 4 + 1) t/|(z)) = 0, 

32q 2 r l+ (z)z 6 +(qz 4 + l) ((qz 4 + l) ( Z V+4) $ 2 (z) - 
(Sg^z)* 2 + (qz 4 + 1) + z&i(z)))) = 0, 

(zV + 4) $ 2 (z) - ^ ($' 2 (z) + = 0, 

(zV + 4) h xy (z) - z (h' xy (z) + zh%(z)) = 0, 

— 32q 2 r] + (z)z 7 + (gwz 4 + u) 2 rj_(z) z— 
- (qz 4 + 1) (8q$ 2 (z)z 5 + (Aq^' 2 (z)z^ + (qz 4 + 1) t£(z)) z + (5g^ 4 - 3) rf_(z)) = 0. 

(119) 

Solutions for the EOM in the Liu-Tseytlin case at zero frequency u = are: 















$4 









/ C 2 (^ 4 + l) 2 + C 1 (^ 4 + l) 



-g 2 C 2 ^ 8 +C3^ 4 +C4 
53 

Ci z 4 +C 9 



V 



C 6 9 2 + (g(^ 4 +2)2 4 +2)(4 (? (Ci+C2)+2C3) 
4 g (<?z 4 +l) 



(120) 



Solution modes for a non-zero frequency: 



28 



$1 = ^K^zu^z 6 + ^^(z^dz 2 , 



§ 2 = C 1 



~iqu) 8 z 10 _ 161quj 8 z 10 
552960 



6144 552960 1 

, go: 8 log(z)z 1Q , qw 8 log^z 10 _ qu 8 log(16)z 10 
H I 6144 92160 



6144 ^ 6144 92160 
qui 8 log(8)z 10 _ qui 8 log(4)z 1Q 1 6 8 _ 169gcj 

27648 184320 192 ' y 2304 



27648 184320 " r 192 ' y 23040 ~ l ~ 

^g^log^)^ 8 + ^g^log^)^ 8 - g^gu; 6 log(16)2 8 - 

gqj 6 log(4)g 8 , 14 6 17 „, ,4,6 i 1 ,4WC ? W 6 4- 

1920 r- jQiqoj z — T^qoJ z + j^qu log{z)z + 

■±qoj A \og(u)z 6 - ^qu 4 \og(4)z 6 + \qu 2 z A ] + \u 2 K 2 (zu)C 2 



$ 3 = -ui 2 K 2 (zu)C 7 



$ 4 = -^^(zu/JC'g, 



1 2 6 1 2 6 4 Zqhiz^dz 2 u 2 K 2 (zu)C 1 z 2 

4> 5 = qui C\z H — qui CaZ — qC\z ; — ; r — - h 

5 12* 6 y y (<^ 4 + l)w 2 qz 4 + l 

4q 2 I 2 (zoo)C 6 z 2 qoo 2 K 2 (zuj)C 6 z 2 
+ (qz 4 + l)co 2 + 8 (gz 4 + 1) 

(121) 

The thermal version of the Liu-Tseytlin backgrounds leads to the following equations 
of motion: 



' (u(u 3 +6u+4w 2 +4fc 2 (?i 2 -l))-3)j7+ + „ 

4«2(n2-l) 2 = U ' 

-4q {u 2 + 1) h{u)u 2 + 4 (u 2 - 1) (2gu/i' + tt 2 {u 2 - 1) 77-") u 2 + 
< W{u{u z + Qu + Au 2 + Ak 2 {u 2 - 1)) -3)77- = 0, (122) 

4 (V (u 2 - l) 2 + 2vr 2 g (2« (u 2 - 1) V + ' - (u 2 + 1) r)+(u))) u 2 + 
k + (« (u 3 + 6u + 4w 2 + 4fc 2 (u 2 - 1)) - 3) h = 0. 
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